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We study the domain number and size distributions in the one-dimensional Ising and g-state 
Potts models subject to zero-temperature Glauber dynamics. The survival probability of a domain, 



S(t) 



t~ 



, and an unreacted domain, Qi (t) ~ t , are characterized by two independent nonrtrivial 



exponents. For the Ising case, we find ip = 0.126 and S = 1.27 using numerical simulations. We 
develop an independent interval approximation (IIA) that predicts the qualitative behavior of the 
domain distribution and provides good estimates for the exponents. Exact results for the domain 
distribution are also obtained in several solvable cases. 
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I. INTRODUCTION 

When a system is quenched from a homogeneous high- 
temperature disordered state into a low-temperature 
state it does not order instantaneously; instead, domains 
of equilibrium ordered phases form on larger and larger 
scales 0. It has been generally confirmed that a scale- 
invariant morphology is developed, i.e. the network of do- 
mains is (statistically) independent of time when lengths 
are rescaled by a single characteristic length scale that 
typically grows algebraically with time. However, even 
for simple coarsening processes little is known about 
more subtle properties such as the domain size distri- 
bution |^,|3|]. One such feature that has attracted con- 
siderable interest recently concerns the "persistence" of 
the local order parameter, the probability that it has not 
changed sign in a given time interval. Persistence has 
been investigated theoretically |^-^| and experimentally 
JTof in spin systems, interacting particles systems [jy] 14 1. 
Lotka-Volterra models [[15 16|, breath figures growth [17|, 
foams , and even simple diffusion |ll],|2(| . 

Single spin persistence provides a natural counterpart 
to the survival probability in the realm of many particle 
systems. In the context of reaction processes, persistence 
is equivalent to the survival of immobile impurities and 
therefore does not provide information about collective 
properties of the bulk. In contrast, domains are the nat- 
ural spatial elements of a coarsening process. In this 
paper, we ask for example what is the survival of an en- 
tire domain, S(t)l This quantity decreases with time as 
a power-law S(t) <~ t~^. Similar to other critical expo- 
nents, ip is universal in the sense that it is independent 
of many details such as the initial conditions. However, 
it is model dependent and in this respect differs from the 
growth exponent that depends only on the conservative 
nature of the dynamics. 

In the present work, we examine systems with short- 
range interactions described by a scalar non-conserved 
order parameter. We focus on the ID Ising model and 
its generalization to the g-state Potts mode l pl| evolving 
according to Glauber spin-flip dynamics |22||. Sequen- 
tial dynamics has been chosen without loss of generality 
as parallel dynamics exhibit similar asymptotic behavior 



p3| . In one-dimensional systems with short-range inter- 
actions, the order-disorder transition takes place at T = 
p4j , and since we are interested in coarsening, we restrict 
attention to zero temperature. 

This paper is organized as follows. In the next sec- 
tion, we discuss the Ising and Potts models. We de- 
fine the domain number distribution and determine it 
analytically in the limiting cases, q — > oo and q — > 1. 
Otherwise, we develop an Independent Interval Approx- 
imation (IIA) that assumes no correlations between ad- 
jacent domains. The IIA predictions compare well with 
Monte Carlo simulations by giving correct description of 
the domain statistics as well as good estimates for the 
underlying exponents. Section II is concluded with the 
q = oo Potts model in arbitrary dimension d > 2. In Sec. 
Ill, we obtain the domain distribution in two solvable 
cases: the Potts model with only energy lowering tran- 
sitions and the deterministic ballistic annihilation model 
16 1 . Section IV discusses some open issues and 



contains a summary. 



II. ISING AND POTTS MODELS 

A. The Ising model and the g-state Potts model 

We start with the ID Ising model subject to T = 
Glauber dynamics p2|] . To examine the role of the num- 
ber of equilibrium phases we also consider a generaliza- 
tion of the Ising model, the g-state Potts model. Experi- 
mental realizations are known for q = 2 (the Ising model) 
and q — 3,4,oo plj. The q = oo case describes several 



polycrystals 



foams 



cellular structures |26[ 
soap froth pi] , and magnetic bubbles j2£f 

We consider uncorrelated initial conditions where each 
of the q phases is present with equal density 1/q. The 
T = Glauber-Potts dynamics proceeds by selecting a 
spin at random and changing its value to that of one of 
its randomly selected neighbors. Thus, domain walls per- 
form a random walk and upon contact, they annihilate or 
coalesce, depending on the state of the corresponding do- 
mains fl3Q-p3]. Identifying a domain wall with a particle, 



1 



(A), and absence of a domain wall with a hole (0), one 
finds the single-species diffusion-reaction process 30 3^] 



.40 



OA 



AA ±$ 00, 



AA ^4 AO or OA. 



(1) 



The rates indicate the relative probabilities by which each 
event occurs. 



B. Domain Number Distribution: Definition and 
Scaling Properties 

Our goal in this study is to investigate S(t), the prob- 
ability that a domain, initially present at the system at 
time t = 0, has not flipped up to time t (see Fig. 1). We 
will present theoretical and numerical evidence support- 
ing an algebraic long time decay of this survival proba- 
bility, 



S{t) ~ t 



(2) 



Such a behavior is robust, as the exponent %p is not sen- 
sitive to initial state of the system (provided long ranged 
correlations are absent). Our results will also strongly 
suggest that the exponent ip is nontrivial, i.e., it can- 
not be extracted from so-far known exponents associated 
with the Ising model. 

+ - + - + + 
T 



m=3 / m=l^ 

Fig. 1. Domain motion in the Ising model. Surviving do- 
mains are marked by +, annihilated domains by — . The do- 
main number at a later time is also indicated. 

In principle, a surviving domain may undergo coales- 
cence with other similar phase domains. Thus, a natu- 
ral generalization of the domain survival probability is 
Qm(t), the density of domains composed of m original 
domains (see Fig. 1). This quantity satisfies the initial 
condition Q m (0) = S m< i. The total domain density, N(t), 
is given by N(t) = J2 m Qm(t), while the domain survival 
probability counts initial domains that have not shrunk 
and hence contains the density Q m {t) with weight m 



s , (*)=E i 



»(*). 



(3) 



The average number of domains contained within a sur- 
viving domain (m(t)) = S(t)/N(t) grows algebraically 



according to (ra(t)) ~ t v ~^ with v the domain decay ex- 
ponent, N(t) ~ t~ v . If the behavior of Q m (t) is truly 
self-similar, it should follow the scaling form 



(4) 



The scaling function Q(z) exhibits the following extremal 
behavior 



z a z«l, 
exp(— kz) z 1. 



(5) 



The small argument tail describes domains that contain 
a very small number of initial domains. In particular, the 
quantity 



i(t)~r 



(6) 



is of special interest: It gives the density of domains 
which avoided merging with their neighboring domains 
up to time t. 

The inequalities Qi{t) < £) m Q m (i) < Y, m m Qm{t) 
lead to the bounds 



tp <v <S. 



(7) 



Taking into account that at least one surviving domain 
surrounds a persistent spin gives P(t) < S(t), where 
P{t) ~ t~ e is the density of persistent spins. Thus we 
arrive at another upper bound 



(8) 



for the exponent ip. We shall show below that these 
bounds are strict for the Potts model and that there are 
models with tp = 9, ip = u, and v = 5. The bounds (0) 
and (^) suggest that a domain decays with the slowest 
rate in the problem. 

A useful relation between the scaling exponents can be 
obtained by substituting m = 1 in Eq. (j|) 



5- v = {v - r/0(l + a )- 



(9) 



Thus, among the three exponents ip, 5, and tr, only two 
are independent. For the g-state Glauber-Potts model 
both the domain decay exponent v — 1/2 and the 

persistence exponent 9(q) H are known. It gives hope 
that analytical determination of the domain exponents is 
also possible. 

Quite obviously, domains disappear when their size 
vanishes, and therefore the domain size and number 
distributions are intimately related. One is therefore 
forced to consider the distribution of domains of size 
n consisting of m original domains at time t, denoted 
by P n ,m{t)- The aforementioned number distribution is 
Qm{t) = J2 n Pn,m{t), and consequently, the domain sur- 
vival probability is S(t) = S„ im mP„, m (t). 

As will be seen later, studying the joint size-number 
distribution requires detailed knowledge of the domain 



2 



size distribution P n (t) = ^ m P n ,m{t). This distribution 
obeys the normalization condition 



where /„ is the modified Bessel function of order n 1341 



,nP n {t) = 1 



(10) 



reflecting the conservation of the total length. The to- 
tal domain density is simply N{t) = ^2 n P n (t)- Since 
the average domain length grows as n ~ t u , the length 
distribution follows the scaling form 



P n {t) ~ t- lv V(nt~ v ) 



(11) 



The length distribution scaling function has the following 
limiting behavior 



V(x) 



X X<1, 

exp(— Ax) x 3> 1. 



(12) 



In subsection D, we shall develop an approximation 
scheme that helps elucidate many of the qualitative and 
quantitative features of the domain size and number dis- 
tributions. 



C. Solvable Cases 

1. The q — > oo Limit 

In the q = oo case, similar phase domains never coa- 
lesce and therefore the domain number is trivial, m = 1. 
Thus P n>m {t) = P n (t)8 m , u N(t) = S{t) = Qt(t), and 
v = ip = 5 = 1/2, The value of the exponent has been 
obtained by noting that domain boundaries perform in- 
dependent random walks and a domain disappears when 
its boundaries meet. Thus domains survive with proba- 
bility identical to that of a random walk in the vicinity of 
a trap, N(t) ~ t~ 1/2 Q, or v = 1/2. On the other hand, 
an individual up spin inside this domain has not changed 
its sign if it has not been crossed by both domain walls. 
Therefore, the persistence probability is proportional to 
t- 1 ' 2 x £-!/2 = t -\ i.e. 9 = 1. Hence, the bound (g) is 
strict. 

The domain length distribution P n (t) obeys the diffu- 
sion equation 



dPn 

dt 



Pn+l + Pn-1 — 2P„, 



(13) 



with the boundary condition Po(t) — 0. This rate equa- 
tion satisfies the length conservation of Eq. ([n]). Solv- 
ing ( Jl3| ) subject to the appropriate initial conditions, 
P n (0) =S n ,i, gives 



P n (t) = [J n _i(2t) - /„+i(2i)] exp(-2i), 



(14) 



and 



N(t) = S(t) = Q 1 (t) = [I„(2t) + h(2t)} exp(-2t), (15) 



The length distribution scales according to Eq. (11), with 
v = 1/2 and V(x) ~ x exp(— x 2 /4) l\pK. The generic ex- 
ponential behavior of Eq. ( |12|) is now replaced by a Gaus- 
sian one, indicating that A — * as q — ► oo. In the long 
time limit, S(t) = Qi(t) = N(t) ~ (Trf)^ 1 / 2 confirming 
the previous heuristic findings v = %jj = 5 = 1/2. 



2. The q — > 1 Limit 

The ID T = Glauber-Potts model with arbitrary 
q > 1 can be mapped to the Ising model with magnetiza- 
tion /j, = 2/q — 1 [p5|| . In other words, the volume fraction 
of the down phase is ip = 1 — 1/q 0. In particular, the 
limit tp — > allows treatment of the limiting case q — ► 1 
by focusing on the majority domains. The typical ini- 
tial size of such domains is ip^ 1 — > oo. This shows that 
in the limiting case q = 1 the minority domains cannot 
meet and the majority domains' sizes change appreciably 
due to coalescence. Thus, majority domains never disap- 
pear, i.e., S(t) — 1 and -0 = 0. Similarly, the persistence 
exponent is found: 9 = 0. A majority domain remains 
unreacted till time t if both of its neighboring minority 
domains survive, Qi(t) = N 2 (t). The density is given by 
the q = oo solution (|l5|), and we find Q±(t) ~ (7rt) _1 and 
5 = 1. 

The number distribution of the majority domains can 
be determined as well. The dynamics proceeds by mi- 
nority domains shrinking to zero and thus leading to co- 
alescence of surrounding majority domains. Such aggre- 
gation events occur independently with rate Pi /N 2 and 
the domain number distribution evolves according to 



dQ„ 



P_ 

N 2 



,i=i 



QjQm-j - 2NQ 



(16) 



subject to the initial conditions Q m (0) — 8 m i- It is 
helpful to absorb the time-dependent rate Pi/N 2 into 
the time variable 



T = f'dt'SMX =N- 1 (t)-l ) 



II 



N 2 (t>) 



(17) 



with the overall density of (|lj|), and the last equality eval- 
uated using N = _p x . With this time variable Eq. @ 
reduces to the classical Smoluchowski equation [ p6| 



~dT 



V QjQm-3 -2NQ ri 



3 = 1 



(18) 



Solving Eq. Qlq ) with the appropriate monodisperse ini- 
tial conditions gives Q m (T) = T"''^ + T)~" 1-1 ~ 
T- 2 exp(-m/T) @. Indeed Q X (T) = N 2 {T) = (1 + 
T)~ 2 , in agreement with the previous argument. In 
the long-time limit, T ~ iV _1 ~ \fwi and Q m (t) ~ 



3 



(7rt) _1 exp [— m(7rt) _1 / 2 ] . Thus the domain number dis- 
tribution scales according to Eq. (^) with the purely ex- 
ponential scaling function 



Q(z) — 7r exp(— ZTT 



(19) 



The average domain properties are S(t) — 1, N(t) ~ 
(nty 1 / 2 , and Q x (t) = N 2 (t) ~ (nt)' 1 and the scaling 
exponents a — ip = 0, v = 1/2, 8 = 1. 

Changes in the domain size due to domain wall diffu- 
sion are negligible here and the joint size-number distri- 
bution evolves according to 



dP n 



dT 



Vp P 

/ j 1 i, J- 1 n— 2,m— j 



2NP n 



(20) 



Eq. (|2^) generalizes the Smoluchowski equations for ag- 
gregation with two conservation laws fl37|] . Introducing 
the generating function F(u, v, T) = ^ n m u n v m P n . m (T) 
one can solve Eq. (20) for arbitrary initial conditions 
to find §7| F(u,v,T) = F (u,v)(l + T)- J [l + T - 
TFq(u, V)\ . In the present case, the appropriate ini- 
tial conditions are P n ,m(0) — 8 m .iip 2 (l — (fi)"^ 1 and 
hence Fq(u, v) = uvtp 2 [l — u(l — tp)]^ 1 . Evaluating the 
limits rup — > n, and t — > oo, we arrive at the scaling 
form P n .m{t) — t~ 5 / 4 <fr(x,y), with the scaling variables 
x = (m + n)(7rt) -1 / 2 , y — (m— n)(7rf) -1 / 4 , and the scal- 
ing function 



$(x, y) = (nx) 1/2 exp(-3j - y 2 /2x). 



(21) 



Instead of the naive scaling variables nt~ x l 2 and mi -1 / 2 , 
unusual scaling variables underlie the scaling function 
(pl|). The former scaling variable x is just the sum of the 
naive scaling variables while the latter "diffusive" scale 
y is hidden. In this case, the domain length and the do- 
main number are equivalent, and their underlying scaling 
functions are identical V{x) = 7r _1 exp(— xn^ 1 / 2 ). Nev- 
ertheless, there is a considerable difference between the 
variables n and m as the latter is generally not a con- 
served quantity. Furthermore, as q decreases from oo to 
1, the decay coefficient A governing the domain length 
distribution increases from to -k^ 1 ! 2 . 

The above results apply for q close to unity as long 
as neighboring domains do not interact, i.e., as long as 
the diffusion time scale is smaller that the domain size, 



9 



Eventually, this no longer holds, and corre- 



lations between majority domains develop. Nevertheless, 
in the limit q — > 1 the equations ([l8]) and ( pp| ) are ex- 
act since no correlations develop if none are present ini- 
tially. Similar reasoning applies to several models where 
domains are immobile and merging occurs 0. 



D. Independent Interval Approximation (IIA) 

Ignoring correlations between neighboring domains al- 
lows us to develop an approximate theory for the time- 
evolution of the domain distribution. This so-called IIA 



proved useful in studies of related reaction-diffusion pro- 
cesses 



1. The Length Distribution 

The joint number distribution requires knowledge of 
the length distribution and we start by deriving a master 
equation for P n (t). Consider first the Ising case. In an 
infinitesimal time interval At, the domain P n {t) changes 
according to 



P»(t + Ai) 

+ MP n - 1 (t) 



(l-2At)P„(f) 
Pi(t) 



1 - 



AtP n {t) 



N(t) 



N(t) 
AiPi(i) 



AtP n+1 (t) (22) 



' 5- N(t) N(t) 



where N(t) = J2 n Pn(t) is the total domain density. The 
first term on the right-hand side of Eq. ([22] ) counts for 
the probability that both domain walls do not hop. The 
next two terms describe gain due to diffusion, with the 
prefactor (1 — Pi/N) in the second term to ensure that 
the hopping domain wall does not disappear. The forth 
term represents loss due to disappearance of the smallest 
domain, located on the boundary of our domain, while 
the final term accounts for gain due to domain merger. 

Eq. J22] ) assumes that the sizes of adjacent domains 
are uncorrelated, and thus is mean-field in nature. In 
the limit At —> the difference equations (^2|) turn into 
a system of differential equations: 



dP n 

dt 



— Pn-l + Pn+1 — 2P n 



(23) 



N 2 



^2P i P n _ 1 _ i -N(P n + P n . 1 ] 



Eqs. ( |23| ) apply for n — 1 if we set Pq = 0. One cru- 
cial test is to verify the length conservation of Eq. ([h]) . 
Another test is to sum all equations in (^) to get 
N = —2P±. This is an exact equation, since three do- 
mains disappear and one is born in each annihilation 
event. 

Generally, in the Potts case, the domain size distribu- 
tion evolves according to the rate equation 



dPn 

dt 



Pn—l + Pn+1 — 2P n 



(24) 



+ 



Pi 



{q-l)N 2 



n-2 



^PiPn-l-i-NiPn + Pn-x) 



Indeed, collision of domain walls results in annihilation 
with probability -^j or in coalescence with probability 

■ Only annihilation events affect the domain distribu- 
tion and thus the — ^ prefactor of the annihilation term. 



4 



In the cases q = 2 and q = oo, Eqs. (|23|),([l3D are clearly 
recovered. In the limit q — > 1, only the reaction term sur- 
vives, in agreement with Eq. (EG). One can also verify 
that the total length is conserved and the total domain 
density decays according to the exact rate equation 



dN 
~dt 



q 



-Pi. 



(25) 



The diffusion term in Eq. ( pij ) implies (n(t)} ~ t 1 / 2 , 
and since (n) ~ A^ _1 the correct decay exponent v = 1/2 
p2| , |30| is recovered. In the following, we will need to 
determine the asymptotic prefactor A : N(t) ~ Ai -1 / 2 , 
A = J dxP(x). The density rate equation ( p5| ) implies 

Pi ~ 7"(0)t- 3 / 2 with P'(0) = a=±A. 

A quantitative analysis of Eq. (|24|) may be carried by 
treating the variable n as continuous. The quantity P(x) 
satisfies 



where V = dV/dx and V * P = f* dyP(y)P(x - y). 
The normalized Laplace transform of the scaling func- 
tion P(x), p(s) = A~ x Jp 00 dx e~ sx P(x), obeys 



The value of the prefactor A predicted by the IIA may 
be compared to the exact one, Aixact = (1 — qT V )l\f^ 
p2[ (see Fig. 2). In the extreme cases of q = 1 and q = oo 
the prefactor A is exact. The mismatch is worst for the 
Ising (q = 2) case where A = r(3/4V/T(l/4) = 0.337989 
while A= xact = (47T)- 1 / 2 = 0.28209 §|. 

0.75 





approximation 




exact 







0.5 
1/q 



Fig. 2. The prefactor A of Eq. @ 

^exact = (1 - q'^/y/n. 



the exact value 



ds qs 



q-2\ q 
. P 

qs 



1 



qs 



(27) 



subject to the boundary condition p(Q) = 1. The trans- 
formation vjs) — l — qs 2 — qsj^ lny(s) reduces the Riccati 
equation ( p7|) into the parabolic cylinder equation, 



ds 2 



1 + - - «* 

q 



y = 0. 



(28) 



The solution to (||) reads y(s) = C-D 1/q (-sV2) + 

C+D 1 / q (s\/2), with D 1 / q (x) the parabolic cylinder func- 
tion of order 1/q Q. The large s behavior of p(s), 
p(s) ~ ^l s - 2 ; implies C_ = 0, and we get 



2q 

p(s) 



l-qs 2 ~qs— \nD 1/q (sV2). (29) 



The normalization condition ^2 n nP n (t) — 1 can be re- 
duced to Ap'(Q) = — 1. This allows us to determine the 
constant 



A = 



£[i- & 



2(/J 



(30) 



where T denotes the gamma function. In deriving ( J30| ) 
we have used the properties |3J] 



D c (x) ~ a; c exp(-a; 2 /4)[l + C(x" 



(31) 



and 
D c (0) = 



^1/22^/2 ^l/2 2 (c+l)/2 

r(i/2- e /2) ' ^ (0) = -^H^T- (32) 



The IIA gives correct qualitative results including: (i) 
density decay N(t) ~ i -1 / 2 ; (ii) linear small size dis- 
tribution, P(x) ~ x ' as can ^ e seen consider- 
ing the large s behavior of p(s) ~ 2g s ~ 2 ' ex P°~ 
nential tail, 7 3 (a;) ~ gAAexp(— As) as in Eq. @. The 
tail follows from the behavior of the Laplace transform 
p(s) ~ q\/(s + A) near its pole at negative s = — A, given 
by the first zero of D\/ q (— AVz) = 0. For the Ising case 
one has A = 0.5409. This value should be compared with 
the exact value A = Q{Z 12) / ky/ir = 0.368468 obtained 
by Derrida and Zeitak 0] and the approximate value 
A = 0.35783 obtained by Alemany and ben-Avraham [||. 



One may wonder regarding the value of the IIA. It is 
analogous to that of Ref. || but the quantitative agree- 
ment is worse in our case. The answer is simple: Our 
approach is self-consistent while the approach of Ref. 0] 
is not. Indeed, making use of the assumption that the do- 
main sizes are uncorrelated, Alemany and ben-Avraham 
express the equal-time two-spin correlation function via 
the domain size distribution P n (t). Then they use the 
exact expression for the equal-time two-spin correlation 
function [^2j to determine P n (t). However, would they 
use their key assumption that the domain sizes are un- 
correlated everywhere they would eventually obtain our 
expression for P n (t). In contrast, our approach is self- 
consistent as all our the results are derived within the 
same scheme. Additionally, to determine the more sub- 
tle characteristics to be described below one does not 
enjoy the luxury of known exact analytical results. 
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2. The Joint Distribution 

We are now in a position to tackle the joint size-number 
distribution, P n . m {t), which captures both the spatial 
and "historical" characteristics of the coarsening domain 
mosaic. The corresponding rate equation is a generaliza- 
tion of Eq. (Mj 



dP n 



dt 



= Pn- 



Pn 



+ l,m 



-2P„ 



(33) 



Pi 



(q-l)N 2 



^ ^ Pj,jPn— 1 — i,m— j ^{,Pn.m ~\~ Pi— l,m) 



with the initial condition P n ,m(0) — <5n,i^m.i and the 
boundary condition Po, m (i) = 0. The variable m is 
"mute" in some sense. It appears in a nontrivial way 
only in the convolution term. One should verify that 
this master equation is self-consistent. First, by sum- 
ming over m, we recover Eq. (|24|). Second, it implies 
that the domain survival probability satisfies the exact 
linear equation dS/dt = -^mfi.m. So far, we have 
not succeeded in solving the joint distribution. Never- 
theless, it is still possible to obtain analytically many 
interesting properties of Eq. (j33|), including the scaling 
exponents. 

Let us consider the distribution of domains which have 
not merged with other domains up to t, R n (t) = P Uj i(t). 
For such domains, the convolution term vanishes and 
they evolve according to the linear rate equation 



dR n 
dt 



Rn-l + Rn+1 — 2i?n 



Pi 



(q - 1)N 



(R n + R n -i) 



(34) 



with the initial condition i? n (0) = 8 n ^\ and the boundary 
condition Ro(t) = 0. In the continuum limit we again 
replace R n -i + R n +i — 2R n by d 2 R/dn 2 and i? n + R n -i 
by 2R n to find a diffusion-convection equation for R n (t) . 
The transformation R n — > R n N~ 2 / q reduces this equa- 
tion to the diffusion equation (^) for R n , which is 
solved to yield R n (t) ~ N 2 / q t- ^(nt" 1 / 2 ), with K(x) = 
xexp(— 2; 2 /4)/y / 7t. The large n behavior of P n ,i{t) mim- 
ics the q = oo case in that it exhibits a Gaussian behavior, 
Pn.i(t) ~ exp(— n 2 /At), while the average domain density 
decays exponentially, P n (t) ~ exp(— n/^/t). Large inter- 
vals are more strongly suppressed when they have a small 
number of ancestors and therefor, the n-tail of the joint 
distribution strongly depends upon m. The total density 
of unreacted domains is Qi(t) = J2 n R n ~ t~^~~i , which 
gives the decay exponent 



5 = 



1 



(35) 



Obtaining the second independent exponent ip is more 
involved. The natural approach, i.e., a direct inves- 
tigation of the domain number distribution, Q m , ap- 
pears to be useless, as it requires knowledge of Pi, m and 



hence the entire P n , m - The domain survival probabil- 
ity can be alternatively obtained by considering U n (t) = 
E m ™^,m(')' This quantity obeys 



dU n 

dt 



U n -i + U n +i — 2U n 



(36) 



P^ 



(q - 1)N 2 



n-2 



2 U i P n 



N{U n + E/n_i) 



obtained by summing Eqs. (p3[). We write U n (t) in a 
scaling form U n {t) ~ t^^^Wjit" 1 ^ 2 ). Asymptotically, 
the domain survival probability reads S(t) ~ Bt~^ with 
B = J dxU(x). The scaling distribution satisfies 



U" +\(xU)' + - )U+^-rU*V = 0. (37) 

2 \ qj qA 



The normalized Laplace transform of the scaling function 
U(x), u(s) = B^ 1 J °° dx e~ sx U(x), obeys 



du ( p(s) + # - 1 \ 2ip 

— = 2 h s ) u , u(0) = 1. (38) 

as \ qs J s 



In deriving ( |38| ) we used the relation U'ifi) = Btp, found 
by integration of Eq. $57\), combined with A = J dxVix). 
Substituting the explicit expression ( [29] ) for p(s) into 
Eq. (H), and solving for u(s) yields 



{ 3 ) = 2^D-f q {sV2) drr~ 2 ^D\ lq {r^i). (39) 



This solution is consistent with the anticipated s — > oo 
behavior, u(s) ~ i(js~ 2 . Furthermore, evaluating Eq. (|39| ) 
near the origin gives u(s) = 1 + F(ip)s 2 ^ + Cs + • • •. 
Therefore for u'(s) to be finite near s = 0, we must have 
F{ip) = 0. Evaluating F(tp) gives 



(40) 



0= / drr-^D 1/q (r)D[ /q (r), 



Interestingly the second domain survival exponent, ip, is 
irrational, in contrast with S. 

For completeness, we write the leading extremal be- 
havior of the function U(x): 



U(x) 



x x — > 0; 

a;exp(— Ax) x — > oo. 



(41) 



This behavior can be easily obtained from the extremal 
behavior of the function u(s). When ,s — > oo, u(s) ~ ,s~ 2 , 
while near the pole s — > — A, one finds u(s) ~ (s + A)~ 2 . 



G 



E. Numerical Results 
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To test the IIA predictions, we performed numerical 
simulations on a spin chain of size L = 10 7 . Random 
initial conditions and periodic boundary conditions were 
used. The simulation data represents an average over 10 
different realizations. For the Ising case, we found the 
exponent values ip = 0.126(1) and 8 — 1.27(2) (see Figs. 
3 and 4) . These values should be compared with the IIA 
predictions of ip = 0.136612 and 6=1. 

As was the case for the persistence exponent, 9, the 
domain exponents strongly depend on q. Numerical val- 
ues of the exponents ip and S are summarized in Table 
1 for representative values of q. As q increases, the ap- 
proximation improves and eventually becomes exact for 
the extreme case q = oo. Thus, ip is overestimated by up 
to 10% and S is underestimated by up to 25%. Hence, 
domains of average number m are better approximated 
in comparison with domains with extremely small m. Al- 
though the estimates are not exact they are still useful 
as they exhibit the correct g-dependence. In the q — ► oo 
limit, the exponents approach the limiting value 1/2 ac- 
cording to 6 = \ + \ derived in Eq. (||j) , and ip = \ — | . 
The leading behavior for ip follows directly from the scal- 
ing relation (^), S and a (oo) = 0. 

In the q — > 1 limit analysis of Eq. (40) suggests that 
ip vanishes according to ip oc (q — 1) . This behavior 
agrees with the q = 1 exact solution and is consistent 
with simulation results of an Ising chain with magneti- 
zation fj, = 2/q— 1. It is practically impossible to obtain 
5 conclusively because unrcactcd domains decay quickly. 
The limiting value as q — > 1 appears to be larger than 
the value suggested by the IIA, 5 = 3/2 — (q — 1). 

We performed several checks to verify that the asymp- 
totic behavior of Eq. (^) and Eq. (H) is robust. For ex- 
ample, it is independent of the initial domain wall con- 
centration (provided that the correlations in the initial 
condition are short range) . We conclude that ip and 6 are 
nontrivial exponent, i.e., they cannot be extracted from 
the known exponents associated with the Ising-Glauber 
model. Similar to the persistence exponent, 9(q), the ex- 
ponents appear to be irrational except for the limiting 
cases q — oo (ip = 8 = h, cr = 0) and, maybe, for q = 2 

(i>=l,S=l<T = l). 

The numerical simulations also confirm that the dis- 
tribution function Q m (t) scales according to Eq. ^ (see 
Fig. 5). The scaling function Q{z), defined in Eq. (||), 
decays exponentially for large argument and is algebraic 
for small argument. The scaling relations combined with 
the simulation values give a = 1.05(5). This is consistent 
with the linear behavior seen in Fig. 5 for z<l. Hence, 
similar scaling functions underlie the domain number and 
size distributions As q increases from 2 to oo, the 

exponent a decreases from 1 to 0, respectively. 
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t 

Fig. 3. The domain survival probability S(t)/S(t — 1) in the 
Potts model. Shown are Monte Carlo simulation results for 
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Fig. 4. The density of unreacted domains Qi(t)/Qi(t — 1) 
in the Potts model. The same Monte Carlo simulation as in 
Fig. 3. 





MC 


Eq. (33) 


q 


ip 


S 


(T 


iP 


S 


2 


0.126 


1.27 


1.05 


0.136612 


1 


3 


0.213 


0.98 


0.67 


0.231139 


5/6 


4 


0.267 


0.85 


0.50 


0.287602 


3/4 


8 


0.367 


0.665 


0.24 


0.385019 


5/8 


50 


0.476 


0.525 


0.03 


0.480274 


13/25 


oo 


1/2 


1/2 





1/2 


1/2 



Table 1: Domain exponents for the g-state Potts model in 
one dimension. Local slopes analysis was applied to the simu- 
lation data. The theoretical ip is from Eq. ( po| ) and 5 = | + ~ . 



7 



o 



0.0 




0.0 1.0 2.0 3.0 4.0 5.0 6.0 



Fig. 5. The scaling distribution Q(z) vs. z for three differ- 
ent times t = 10 2 , 10 3 , 10 4 . The data represents an average 
over 100 systems of size L — 10 5 . The inset demonstrates the 
exponential behavior of the large- z tail. 

Direct numerical integration of Eq. (^3|) reveals a num- 
ber distribution, Q m (t), that scales according to Eq. (||), 
and has an exponential tail in agreement with the simula- 
tion results. Moreover, the emerging S(t) falls within 5% 
of the actual survival probability over a significant tem- 
poral range, t < 10 3 . In summary, in addition to predict- 
ing the correct scaling behavior, Eq. ( pS) ) provides a good 
approximation for many quantitative features of the do- 
main distribution, and in particular, good estimates for 
the decay exponents. 



spin in the sea of down spins. On the other hand, for the 
Potts model with symmetric initial conditions the den- 
sity of any phase is globally conserved. This suggests a 
correspondence between the Potts model and the Ising 
model with globally conserved dynamics. It appears dif- 
ficult to make such correspondence rigorous, although it 
is supported by several tests [ pSf , The reduction to the 
Ising model with magnetization fj, = 2/q — 1 in d dimen- 
sions can be hardly considered as a simplification except 
for the q — oo case which may be analyzed within the 
framework of the Lifshitz-Slyozov theory jij. Indeed, as 
q — ► oo the minority (up) phase approaches infinitesi- 
mal concentration, so up domains do not interact and 
the Lifshitz- Slyozov approach, suitably modified for the 
present case |35| . should be exact. 

So consider a set of bubbles (the domains tend to be- 
come round in high dimensions) of up phase in a sea of 
down phase. We cannot restrict ourselves to the single- 
domain situation as in one dimension since when d > 1 
a single bubble would not evolve, while the set of bub- 
bles do evolve: small bubbles shrink and large bubbles 
grow. We are not interested in details of the bubble 
evolution; the only relevant feature is that the radii 
distribution scales asymptotically, N(R,t) N(R,t) ~ 
n-( d+l W(R/1l), with the average radius U ~ \fl. This 
behavior is due to the nature of non-conserved dynamics 
|jj and the prefactor guarantees a conserved magnetiza- 
tion. Clearly, S(t) is determined by computing the num- 
ber of surviving bubbles, S(t) ~ / dRN(R,t) ~ K- d , 
implying 

%jj = 5 = e = d/2. (42) 



F. q = oo Potts Model in Higher Dimensions 

So far, our discussion has been restricted to one dimen- 
sion. Domains are not necessarily well-defined in higher 
dimensions. In particular, in the Ising case it is not clear 
whether our results can be properly generalized to higher 
dimensions. However, in the q — > oo limit of the Potts 
model domains are well defined and dynamics consid- 
erably simplifies |35| . In this limit, it has been argued 
heuristically and confirmed numerically J38|,[l8| and ex- 
perimentally (3^] for evolving soap froth, that 8 = 1 in 
two dimensions. 

We now present a simple heuristic argument which 
gives the exponents tj) = 8 = 9 = d/2 for the q = 00 
Potts model in d > 2 dimensions. First we note that 
an exact correspondence between the dynamics of the 
q-state Potts model and the Ising model with magneti- 
zation fi = 2/q — 1 holds only in one dimension. This 
happens due to the global conservation of the magnetiza- 
tion |^] , a peculiar property of the ID Ising model with 
zero-temperature Glauber dynamics. This global conser- 
vation for the locally non-conserved dynamics does not 
happen when d > 1, as seen by considering a single up 



III. EXACTLY SOLVABLE CASES 

Given that obtaining the exact domain distribution in 
the g-state Potts model appears to be a difficult problem, 
it might prove useful to study simpler problems which are 
exactly solvable. We present in this section exact results 
for a variant of the Potts model with simplified dynamics 
and for ballistic annihilation. 



A. Diffusionless Dynamics 

Consider the T = q-state Potts with simplified dy- 
namics where only energy lowering transitions are al- 
lowed. Thus, domain wall diffusion AO # OA in Eq. (|l|) 
is forbidden and the reaction scheme is 

AA ^4 00, (43) 

AA ^\ AO or OA. 

When q — 2, exchanging the roles of domain walls and 
vacant sites, this problem is equivalent to random se- 
quential adsorption of dimers. Similarly, the q — 00 case 
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reduces to monomer adsorption subject to a volume ex- 
clusion constrain |4C|] . 

Assuming that neighboring intervals are uncorrelated, 
the domain length density rate equation reads 



dP» _ Pi 
dt (q - 1)N 2 



2NP n 



(g - 2)Pi 
(q-l)N 



[Pn-l ~~ Pn] ~ S nt iPi. 



(44) 



For simplicity, we consider the antiferromagnetic initial 
condition P n (Q) — 6 n ,i- While the annihilation term 
is similar to Eq. (pi]), coalescence events are no longer 
offset by domain wall diffusion and thus the second 
term which is proportional to the coalescence probabil- 
ity (q — 2)/(q — 1). It can verified that the total length 
J2 n n Pn = 1 is conserved. 

The density decays according to the familiar rate equa- 
tion ( |25| ) N = — zj-Pi- O n the other hand, the minimal 



gap density satisfies P\ = —P\ 



1 



qPi 



It is useful 



(q-l)N\ 

to introduce the normalized quantity p\ = P\/N which 
obeys p\ = —p\ and thus pi = e -t . Using this result, the 
minimal gap density and hence the density is found 



tf(f)=exp[-«(l- e -*)/(g-l)]. 



(45) 



This agrees with known exact results for the q = 2 
and q = oo cases Usually in random sequen- 

tial adsorption problems, it is convenient to study the 
complementary density of gaps between domains which 
satisfies a linear rate equation. Nevertheless, the IIA 
is exact in this case as no "mixing" of domains due to 
diffusion occurs. The final domain density is given by 
N(oo) = exp[— q/(q— 1)]. The systems quickly reaches a 
jamming configuration where domain walls are isolated 
and immobile. Thus, no coarsening occurs and the pos- 
tulated scaling behavior does not apply. Nevertheless, as 
will be shown below, the domain size and number dis- 
tributions and in particular their tails do resemble their 
diffusive counterparts. 

The length distribution P n can be found using normal- 
ized distribution p n = P n /N which satisfies (for n > 2) 



dpn 
dt 



Pi 
q-1 



n-2 

E 



PiPn-l-i + {q - 2)p n - 



(46) 



This equation can be further simplified by introducing 
the modified time variable T, defined via dT/dt = p\ im- 
plying T = 1 — e _t . To solve Eq. (^) we introduce the 
generating functions 



(47) 



that satisfies 



*|^ = -^[ P (z,r) 2 + ( 9 -2M z ,r)]-z. (48) 



Solving Eq. (E§) subject to the monodisperse initial con- 
ditions P n (0) — S n ,i, i.e., p(z,T = 0) = z, we get 



p(z,T) = l+ q 



z + q-1 



exp 





- 1 







(49) 



Clearly, quantities such as the domain density and the 
domain length distribution approach exponentially fast 
their limiting values. We are especially interested in these 
limiting values. By expanding the generating functions in 
powers of z we get the limiting density of short domains 

P»M = °' 2^J» and i0TH = 1 ' 2 ' 3 ' 4 ' 

Similar to the behavior seen in the Glauber-Potts 
model, large domains are suppressed exponentially, 



P„(t) ~[A(g,f)] f 



n > 1. 



(50) 



Here \{q,t) is equal to the inverse of the first simple 
pole of the generating functions p{z,t). This can be 
seen from the ^2i n X n z n oc (A -1 — z)^ 1 . Note that 
A(g, i) vanishes when t — > as implied by the ini- 
tial conditions. When t — * oo, X(q,t) — > Xoo(o)- In 
particular when q — > oo, A 00 ((3 , ) vanishes according to 
Aoo (q) ~ 1/ (In q) , indicating a faster than exponential 
limiting behavior. Indeed, when q = oo the generating 
functions p(z, oo) = l + (z— l)e z gives an inverse factorial 



decay, p n (oo) = (n- l)/n! . 



-n In n 



. In the complemen- 



tary g=l limit, the density vanishes, as was the case in 
the diffusive counterpart. 

The domain number-size distribution can be obtained 
by generalizing Eq. (E3) 



dP„ 



Pi 



dt (q - l)N 2 
(g-2)Pi 



^ ' y , Pj,j Pn—X—i,m—j 2NP n 7n 



[Pn-l,m — Pn,m] — ^n,l^m,\Pl- (51) 



A solution using the generating functions technique is 
possible here as well. However, this solution is too cum- 
bersome, and we briefly discuss its qualitative features. 
There are two limiting cases. When q = oo, the joint 
domain- number distribution simplifies to P nm — P n 5m,i- 
When q = 2, domains are always of odd length and 
Pn,m — Pnfim,(n+V) /2- Hence, the domain number dis- 
tribution also decays exponentially, Q m {t) ~ [A(q,t)] n . 
Similar to the length distribution, the decay constant 
A(q) vanishes when q — > oo. 

In summary, although the restricted dynamics Potts 
model does not exhibit coarsening or scale invariance, 
the number and size distributions mimic their diffusive 
counterpart large n and large q behavior. 
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B. Ballistic Annihilation Model 

Consider a binary reaction process with particles mov- 
ing ballistically and annihilating upon collision. Assum- 
ing a bimodal velocity distribution, we set these velocities 
equal to ±1, without loss of generality. Identifying do- 
main walls with particles, this two- velocity ballistic anni- 
hilation process p5[ is equivalent to deterministic coars- 
ening in a system with three equilibrium states [p"5| , [T6[ . 

The domain size distribution for this ballistic anni- 
hilation process has been investigated in f[q| . Here 
we want to compute the domain survival probability 
S(t). There are actually four such survival probabili- 
ties depending on the initial velocities of boundary in- 
terfaces; we denote the corresponding survival proba- 
bilities by £+_(*), S_+(t), and S—(t). Then 
the total survival probability is just the sum S(t) = 
±[S ++ (t)+S + -(t) + S-+(t)+S—(t)}. We need to 
specify the initial conditions. Let us assume that in- 
terfaces are located according to the Poisson distribution 
with unit density. For such symmetric initial conditions 
we have S++(t) = S--(t). One immediately gets the sur- 
vival probability in the simplest case when the interfaces 
move toward each other: 



procedure so that Wj 



Wi. 



-21 



(52) 



To compute the survival probability of parallel mov- 
ing interfaces we note that the probability S+(t) for 
a single right-moving interface to survive is S+(t) = 
[S++(t) + 5-| (t)]/2. Combining the known result E3] 



S+(t) = SL(t) = e 
with Eq. (p2T) we arrive at 



it 



7 (2t)+/i(2f)], 



(53) 



S++(t) = S—(t) = e- 2t [2I (2t) + 2h(2t) - 1]. (54) 

" J»(t) ^ 



Using the asymptotic behavior [ 
when t — > oo, we get S++(t) ~ 2/y/nt. 

We turn now to the more challenging problem, i.e. to 
computation of the survival probability when the inter- 
faces move away from each other. The final answer is 
relatively simple: 

S-+ (t) = e- 2t [2Ji (2t) + Ah (2t) + 2I 3 (2t) + 1] (55) 

so that the total survival probability reads 

S(t) = ^e- 2t [2I (2t) + 3I 1 {2t) + 2I 2 (2t) + I 3 {2t)}. (56) 

Eqs. (H|) and (|5^) imply an algebraic asymptotic behav- 
ior, S-+(t) ~ and S(t) ~ 2/V^t for t -> oo. 
Thus, the domain decay exponent is ip — 1/2. 

To derive Eqs. (55) and ( [Sq ) let us consider a sequence 
of interfaces starting from the right interface of our do- 
main as a random walk. Namely, we set Wq = 0, then 
define W\ = Wo + v\ = Vi where Vi = 1 is the veloc- 
ity of the right interface of our domain. We repeat this 



Vj and we treat Wj as 



the displacement of the random walk, started from the 
origin, at the j th step. When the displacement becomes 
negative for the first time, the corresponding interface 
will move to the left and will eventually destroy the do- 
main. Let we meet this interface after 2 N + 1 steps. The 
corresponding probability Pn is readily determined by 
random walk methods [ES3f: 



P 



N 



%-2N 



(27V)! 



(JV + 1)!JV! 



(57) 



The same analysis applies to the left interface of the do- 
main. Thus we have 2N + 1 interfaces to the right and 
2M + 1 interfaces to the left. Our original domain sur- 
vives till time t if the distance between the extreme in- 
terfaces is greater than 2t. In other words, the interval 
of length 2t with the left boundary at the initial location 
of the extreme left interface should contain 2N + 2M 
interfaces at most. The probability of this event is 



U. 



N+M 



-21 



N+M 

E 

fc=0 



{2t) k 



.2N+2M 



due- u - -. 

(2AT + 2M)! 



(58) 



(59) 



The survival probability S \- (t) is now given by 

5 '-+W = Yl P NPMU N+M (t). 

N,M>1 



It proves convenient to expand the summation in (|59j) to 
N = and M = 0. This gives 



S-+(t)= PvPm^jv+mW 

N,M>0 

-2Y / PNU N {t) + U (t). 



(60) 



N>0 



The second sum in Eq. mu) can be rewritten as 



00 2 AT 
U 



poo 

Y^P N U N {t)= / due- u J2 p N 

N=0 •' 2t N=0 



(2N)\ 



(61) 



The sum on the right-hand side of Eq. (^Tj) is compressed 
into J2n>o P n u2N /( 2N )- = 2Ji(u)/u, and the resulting 
integral is 



due -uIlM. =e- 2t [J (2t)+/ 1 (2t)l. (62) 

2t U 



Thus S \-(t) becomes 

poo 

S-+(t) = / due~ u p N p i 

•' 2t N,M>0 

-4e- 2t [/ (2i) + /i(2<)]+e- 2t . (63) 



.2N+2M 



— M (2A + 2Af)! 

N,M>0 v 1 
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S(t) is found from (§§), (g), and (||) to yield 

.2N+2M 



1 / >0< 



^ PnPm 



N,M>0 



(2N + 2M)\ 



(64) 



To perform the summation in Eq. 
binatorial identity ^2 



N+M=L 



PnPm 



|) we need the com- 
= 4P^ +1 which can 
be checked directly. One can also establish this identity 
geometrically by noting that 2 2N P N = (2N)\/N\(N + 1)! 
gives the number of random walks starting at the origin 
and returning to the origin for the first time after 2N + 2 
steps |33|] . An appropriate counting of all such walks of 
length 2L + 4 then leads to the above identity. Making 
use of this identity we reduce Eq. (B4) to 



se>-ijf>.-E(j) 

due- 11 



u\ 2L 2(2L+1) 
,-„ II ( rT 2)! 
h(u) 3I 2 (u)' 



(65) 



due _u [J (u) - J 4 (u)] 



In deriving the second line we have used the definition 
of the modified Bessel functions; the third line has been 
derived by applying the identity [Q J n _i(tt) — I n +i i u ) = 
^ /„ (u) . Computing now the integral in the last line of 
Eq. (|65|) we arrive at ([56|). This completes the proof of 
Eqs. (||) and ©. 

One can try to compute Q m (t), the domain number 
density. First we note that Qi(t) with specified bound- 
ary velocities can be readily found: Qh (t) = e~ 2t , 

while other single-domain densities Q++(i) = Q__(t) 

and Q can be expressed via single-particle survival 

probabilities S±(t): 



. + (t) = S+(t), Q-+(t) = S-(t)S+(t). 



(66) 



Then the total single-domain survival probability reads 
Qi{t) = [Q ++ (i)+Q + _(i)+Q_+(t)+Q__(t)]/4. Asymp- 
totically, Q\it) ~ (Ant) 1 / 2 implying exponents 5 — ip — 
v = 1/2. Moreover the persistence exponent 9=1, and 
thus all exponents are identical with those of the q = oo 
Potts model. These two models exhibit several other sim- 
ilarities p4| . However, the present deterministic model 
of coarsening is quite different in that the number dis- 
tribution Q m (t) is nontrivial. The determination of this 
distribution is more involved and left for the future. 



IV. DISCUSSION AND SUMMARY 

Even in one dimension there are many interesting sit- 
uations where the above coarsening exponents are un- 
known. The simplest case is a diffusion equation which 
can describe coarsening in systems with non-conserved 
order parameter Q . Recently, the persistence character- 
istics for the diffusion equation process have been investi- 
gated M] numerically and theoretically by an approach 



close to the II A. Given the enormous role played by the 
diffusion equation in science, surprisingly little is known 
about its underlying coarsening process p(J . 

Another well known coarsening process is the ID time- 
dependent Ginzburg-Landau equation for a scalar non- 
conserved order parameter H], In this system, domains 
do not move and the coarsening proceeds via flipping of 
the shortest domains. The minimum size grows loga- 
rithmically flsj ], so it is convenient to define the coars- 
ening exponents in terms of the minimum size L rather 
than time t. This process is solvable in that the domain 
size distribution, P n {L) = L- 2 V(n/L), is known fHfj]|. 
The same expression holds for the domain number dis- 
tribution. Some coarsening exponents are simple, v = 1 
and S = a = oo. In contrast, the persistence exponent, 
8 = 0.1750758 Q, is non-trivial. This process resembles 
the q — > 1 Potts model where ip = 6 as well. 

It would be interesting to extend of our work to coars- 
ening systems with conserved order parameter. Besides 
the dynamical exponent z = 3 little is known even for 
the one-dimensional Ising model with Kawasaki spin- 
exchange dynamics [ff6"f . The coarsening exponents ip 
and (5 appear to be non-trivial for the Ising-Kawasaki 
model J47| as well. Another possible direction is to study 
the coarsening exponents for the natural generalization 
of the ballistic annihilation process, the iV-species Lotka- 
Volterra process |l^,Q . 

In summary, we introduced the domain size distribu- 
tion and showed that it obeys scaling and is character- 
ized by two independent nontrivial decay exponents. The 
survival probabilities of a domain and an unreacted do- 
main are described by the exponents ip and S, respec- 
tively. Generally, these exponents obey < ip < 8 and 
ip < v < S. In the most examples the above inequalities 
are strict; however, there are counter examples where 
tp = 6 and/or 5 = v. For the ID T = g-state Potts- 
Glauber model we developed the IIA that predicts the 
correct qualitative behavior of the domain size and num- 
ber distributions and even reasonable estimates for the 
decay exponents. We also worked out the analytically 
tractable limits of q — > 1 and q — > oo. It still remains, 
however, to obtain the exact behavior for general q. This 
might be possible using the techniques used in studies 
of single-spin persistence [|]||,[lj] • I n a static version of 
the Potts model, an exact solution was presented and 
exponential decay of the domain density still occurred. 
It was also shown analytically that the coarsening expo- 
nents in the solvable deterministic ballistic annihilation 
model and the stochastic q = oo Potts model are identi- 
cal. 

These results indicate that several nontrivial decay 
laws underlie the evolution of elementary processes such 
as the non-equilibrium Ising model. These nontrivial ex- 
ponents do not emerge naturally from studies of tradi- 
tional quantities such as spatiotemporal correlations. It 
remains a challenge to find and obtain these underlying 
"hidden" exponents from a more systematic method. It 
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is also intriguing to determine whether an entire hierar- 
chy or a finite number of independent decay modes are 
present in these systems. 
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